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Spatially explicit models have been widely used in today's mathematical ecology and epidemiology 
to study persistence and extinction of populations as well as their spatial patterns. Here we extend 
the earlier work-static dispersal between neighbouring individuals to mobility of individuals as well 
as multi-patches environment. As is commonly found, the basic reproductive ratio is maximized 
for the evolutionary stable strategy (ESS) on diseases' persistence in mean-field theory. This has 
important implications, as it implies that for a wide range of parameters that infection rate will tend 
maximum. This is opposite with present results obtained in spatial explicit models that infection 
rate is limited by upper bound. We observe the emergence of trade-offs of extinction and persistence 
on the parameters of the infection period and infection rate and show the extinction time having 
a linear relationship with respect to system size. We further find that the higher mobility can 
pronouncedly promote the persistence of spread of epidemics, i.e., the phase transition occurs from 
extinction domain to persistence domain, and the spirals' wavelength increases as the mobility 
increasing and ultimately, it will saturate at a certain value. Furthermore, for multi-patches case, 
we find that the lower coupling strength leads to anti-phase oscillation of infected fraction, while 
higher coupling strength corresponds to in-phase oscillation. 

PACS numbers: 87.23.Cc, 05.40.-a, 82.40.Ck 



I. INTRODUCTION 

In recent years, there has been an increased aware- 
ness of the threats posed by newly emerging and high- 
profile infections disease, such as SARS, the H5N1 strain 
of avian influenza, HIV, Ebola, Whooping cough, Dengue 
fever and the spread of influenza. These diseases exhibit 
large scale spatial contagions and long-term spatial pat- 
terns ;i, 2, 3, 4, 5, 6, 7, 8, 9]. 

At present, there is a great deal of interest in the role of 
spatial structure in both ecology and epidemics. All sys- 
tems are to some degree spatially extended; however, the 
classical theory of the dynamics of epidemics ignores spa- 
tial effects with its assumption of homogeneous mixing. 
It is likely that local processes will play an important 
role in the majority of infectious diseases' interactions 
particularly when an infection occurs through the direct 
contact of infected and susceptible individuals [l3, [IH- 
As a consequence, there has been some attempt to take 
explicitly spatial structure into account when examining 
the spreading of parasites [I, [H [H [H, [H [H, [13 ■ For 
example, recent models show that interactions between 
local populations r nay generate complex spatial patterns 
by dispersal [l|, 0, O, [11, [ll|. In modern societies, in- 
dividuals can easily travel over a wide range of spatial 
scales. The interconnections of areas and populations 
through various means of transport have important ef- 
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fects on the geographical spread of epidemics [2| ■ In par- 
ticular, the spatial structure and the different levels of 
the diffusion and transport processes are responsible for 
the heterogeneity, that is an erratic outbreak patterns ob- 
served in the worldwide propagation and persistence of 
diseases (HjUjlElilillS, recently documented for synchro- 
nization and waves 0, H, d, [H, [2 ^ . [2ll and large-scale 
spatial patterns (e.g. spiral wavesp]. llSl. [2^ and Turing 
patterns [2l,[23|) in measles, dengue fever, SARS, and in- 
fluenza. In order to describe such a complex phenomenon 
and develop powerful numerical forecasting tools, differ- 
ent levels of description are possible, ranging from a sim- 
ple global mean-field to detailed individual-based simu- 
lation (see Refs. [T^ l25l l26l. l27l. I2R [29l. Isol and references 
therein) , cellular automata |ll Il4 Il5ri22l. [sil [3^ . [H, H^l , 
coupled map lattices [SSj, etc. 

Dispersing individuals may react in a complex manner 
to local ecological and epidemiological situations P.[5l[2^. 
m, [131, but rare long-distance dispersal events may be 
important in nature [j, ,22, 38, 39]. The modeling meth- 
ods that summarize dispersal in a diffusion coefficient in 
the classical theory-partial differential equations-cannot 
give insights into the importance of rare events . For 
humans and other social animals in which hosts are dis- 
tributed in heterogeneous patches on a large scale, there 
are two critical sides to transmit. The first is local trans- 
mission among individuals within patches or communi- 
ties (cites, towns, and villages). The second is the trans- 
mission between patches. With this point of view, there 
is a clear conceptual link between ecological and epidemi- 
ological theory; a body of recent works have focused on 
the analogies of the spatial dynamics between infectious 
disease and ecological metapopulations 0, H, HO, Elj . The 
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twofold issues are to dissect how infection processes at 
the local scale determine spatio-temporal patterns of epi- 
demics and understand how these patterns are affected by 
the spatial spread between neighborhoods or patches, for 
instance, the epidemic wavefronts observed in the spatio- 
temporal spread of the Black Death in Europe from 1347 
to 1350 di,!!! and West Nile virus [43|. 

Generally speaking, several generic questions can be 
asked when studying the dynamics of epidemic sprea ding . 
One is the explanation of the possible oscillations [l05j | 
with respect to the temporal evolution of the densities 
of infection waves in homogeneous or heterogeneous re- 
alistic situation, as well as the spatial resonance prob- 
lem. Another one concerns the study of the possibility 
that a disease spreading eventually arrives at a steady 
states and the time arrived [1^, In addition, our 

model can also investigate what enables some species to 
persist while others become extinct. This question has 
shaped the history of research on population dynamics 
and remains a central issue for ecologists and epidemiol- 
ogists 

Although the spatial structure has been developed for 
susceptible-infected-resistant (SIRS) models in previous 
studies, the modulation of such effects in combination 
with individuals' mobility is also unknown. We aim to 
bridge this gap in this paper. Specifically, we aim to elu- 
cidate the phase space on the extinction and persistence 
and contrast them with the prediction obtained mean- 
field, as well as the effect of mobility on persistence crite- 
ria. We also examine how the spatial pattern depend on 
the mobility and coupling strength within multi-patches 
case. 



II. MODEL AND METHODS 

A. Cellular automata model 

In a recent report, van Ballegooijen and Boerlijst 
(vBB) showed that there were three types of spatial pat- 
tern in a spatial susceptible-infected-resistant model for 
disease dynamics by using a grid-structured contact net- 
work, also called cellular automata model jH . They were 
localized disease outbreaks with self-limiting in size, tur- 
bulent waves, and stable spiral waves respectively. Fur- 
thermore, they predicted that there existed a trade-off 
between the parameters of infection period and infec- 
tion rate, which emerges from the evolutionary dynamics 
of the system on the stable spiral waves, and referred 
to this relationship as emergent trade-off. These results 
give a guide to understand the spatial features of epi- 
demic such as wave speed and wavelength, as well as the 
persistence and extinction relationship depending on the 
parameters. To model the spatio-temporal phenomena, a 
feasible approach is the combination of grid-based mod- 
els and movements of individuals. In this paper we take 
the framework proposed by vBB as a starting point in 
assessing the phase transitions between the persistence 



and extinction rather than emergence of the stable spiral 
waves. 

We use a spatial susceptible-infected-resistant model. 
The population with N individuals is categorized accord- 
ing to its infection status: susceptibility (S), infectious 
(I), or resistant (R) [U, i, HO, El . Within a subpopula- 
tion, the dynamics for the local populations obey a ba- 
sic reaction scheme conserving the number of population, 
which has been studied both in physics and mathematical 
epidemiology, namely the stochastic infection dynamics 
process identified by the following set of reactions [l| 

S + I^I + I, (la) 
I ^R, (lb) 

R^S. (Ic) 

The first reaction (Eq. (fla|l ) reflects the fact that an in- 
fected host (I) can infect susceptible (S) neighbors with 
infection rate /3; the second reaction (Eq. (HQ) indi- 
cates the acquisition of resistance that hosts are infec- 
tious for a fixed period r/, after which they become 
resistant (R); and the third reaction (Eq. (fTc| ') indi- 
cates the loss of resistance that after a fixed period 
T/j, resistant hosts once again become susceptible. In 
the sense of individual/local- level model, infected hosts 
can infect adjacent susceptible hosts at infection rate 
f3. Here, the infection rate is not equal the proba- 
bility of infection, see P, HI, |46j. In the short-range 
case, the infectious neighborhood consists of two/eight 
direct nearest-neighbors (NN) in a 1D/2D square lat- 
tices. In addition, some recent investigations have in- 
cluded long-range processes in one- and two-dimensional 
lattice models, reporting quite different results on the 
occurrence of self-sustained oscillation, extinction and 
persistence [13, H IM US HII, [H, [H, [H, HI, as well 
as increasing coupling between subpopulation (subcom- 
munity) in spatially structured environment can lead to 
population outbreaks 0, SH^- Thus, an understanding 
of the joint effect of short- and long-range interaction on 
the individuals (including the mobility called as the diffu- 
sion, hopping, or stirring, hereafter, we refer as mobility) 
is desirable from ecological and statistical physical point 
of view ^57i] . Of course, in present paper, we also further 
study the outbreak of disease inspired stochastic lattice 
epidemic model with a nearest-neighbor interaction. 

For ease of comparison, we follow the vBB's model by 
considering a regular network of sites, each of which con- 
tains one of a single susceptible individual (S), an infected 
individual (I) and resistant (R). The susceptibles are in- 
fected by contact with an infected host at a rate f3, and 
the transmission can only occur locally. Every time-step 
At, cells' state change according to the rules in earlier 
studies 0, Q , which are illustrated in the corresponding 
Eq. C]). 
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B. Rate equations 

The deterministic rate (or mean-field) equations de- 
scribe the temporal evolution of the stochastic lattice sys- 
tem, defined by the reactions ([1]), in a mean- field context, 
i.e. they neglect all spatial correlations. They may be 
seen as a deterministic description (for example emerg- 
ing in the limit of large system sizes) of systems without 
spatial structure. The study of the rate equation is the 
ground on which the analysis of the epidemics spread- 
ing and outbreak. In particular, the properties of the 
rate equations are extremely useful for the epidemics per- 
sistence by estimated values for the classical expression 
of 7?.o-the basic reproduction number. TZq is defined as 
the number of secondary infections caused by a single 
infective during its infectiousness period in an entirely 
susceptible population [H, |4^. For the standard SIRS 
mean- field approximation, also referred in reaction ([T]), 
are as follows: 

^ = -l3SI + TnR, (2a) 



-j^=PSI- Tjl, (2b) 

where S and / denote the number of susceptibles and in- 
fected, respectively. The number of recovered individual 
R is obtained by conservation of the entire population, 
i.e., R{t) ^ N-S{t)-I{t). The key quality describing the 
infection is the basic reproduction number TZo = P/ti. If 
TZq > 1 and the initial relative number of susceptibles are 
greater than a critical value Sc — an epidemic de- 

velops {dl/dt > 0). As the number of infected individu- 
als increases, the number of susceptibles S decreases, and 
thus the number of contacts of infected individuals with 
susceptibles decreases until S — Sc, when the epidemic 
reaches its maximum and subsequently decays. This pro- 
cesses alternative again with the time going if TZq > 1. 
Fadeout is most likely to occur in the period immediately 
following a major epidemic (the so-called inter-epidemic 
trough). An important observation is that deterministic 
models cannot reproduce such effects: the ability of pop- 
ulations to recover from very low levels is a well-known 
weakness of many deterministic models [1, [s^ . 

III. SIMULATION AND RESULTS 

To begin with, we revisit some of the findings by vBB 
[l|. Their results then serve as a reference point for il- 
lustrating the difference of our elaborated approach from 
their results. This is justified by the fact that vBB only 
provide the spatial pattern resulting from static dispersal 
among individuals within single patch, as well as emer- 
gence of trade-offs between spiral waves and turbulent 
waves. Our crucial measures are the phase transition 
between the global extinction and persistence of the epi- 
demic, and the effect from the movements of individuals 



on it. In these simulations we set there is no migra- 
tion among patches, i.e., coupling parameter LR = 0, 
in which our simulations faithfully repeat the general 
findings by VBB fl^ that three types of spatial patterns 
emerge (see Fig.[T]) but find the differences when the mo- 
bility of individuals is turned on. 



A. Phase transitions of persistence and extinction 

In this subsection, we analyze the relationship of in- 
fection rate and infection periodic on phase transition 
between global persistence and extinction within single 
patch. In the final subsection, we will introduce the mul- 
tipatches into the spatial structure model by the coupling 
strength between the two patches. 

When the individuals are statical in the local site, then 
the infection only happens in the nearest neighbors, this 
form referred to as local interaction. Specially, we can 
study the dynamics of such spatially extended model in 
one- or two-dimensional space by defining appropriate 
microscopic rules. Equation ([T]) shows the explicit tran- 
sition rates of cells' state in space, where the susceptibles 
infected by infectious is isotropic. Grid size used in the 
simulations is 200 x 200 cells and time-step equals to 
0.01. Larger grid sizes do not change the evolutionary 
dynamics (see the computer code for the CA model at 
http : / / 7y. nuc . edu . cn /j inzhen/English 1 . aspx) . 

The persistence phase space is that the epidemics will 
persist when the parameters lie in this domain. Here, we 
refer to (t/, /9)-parameters space. Their values are de- 
termined by simulating 200 x 200 cell lattices with ten 
independent runnings. Each lattice is initially set with 
random occupied by 100 infected individuals. For each 
obtained data, we fix the infection period r/, and chose 
an initial infection rate (3. Then, increased (3 by small 
step (we use A/3 = 0.05) and 10 thousand time units 
are simulated with that infection rate value to determine 
the critical value of the phase transition. Here, we find 
that 10 thousand time units is a long-term run for per- 
sistence on this system. However, the typical waiting 
time T{N) until extinction occurs is generally very long 
when the system size N is large. This suggests to con- 
sider the dependence of the waiting time T{N) on N. 
Quantitatively, we discriminate between persistence and 
extinction by using the concept of extensively, adapted 
from statistical physics [5^, [BO] . If the epidemic is not 
extinct at that time, the pair parameters' values are as- 
sumed to be within the persistence space. In the present 
model, we find that the epidemic will be persistence when 
parameter (3 is between a minimal critical infection rate 
(3min and a maximal critical infection rate f3rnax ■ It is also 
worth noticing that the definitions of persistence and ex- 
tinction in the presence of absorbing states are intimately 
related to the concept of long-term. 

The threshold result for the deterministic model can 
be used to identify two domains in the space of essential 
parameters TZq, namely TZq < and TZq > 0, where the 
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model solutions behave in qualitatively different ways. 
We endeavor to describe the spatially stochastic model 
in a similar way, by identifying domains in the space of 
its essential parameters where model solutions behave in 
qualitatively different ways. We perform extensive com- 
puter simulations of the system ([T]) on the phase transi- 
tion, which describes the persistence and extinction with 
respect to the (r/, /3)-parameters space shown in Fig. [2] 
which also demonstrates that there exists a trade-off be- 
tween infection period and infection rate. In order to 
compare the classic mean-field approximation of non- 
spatial theory TZq and the spatial structure of individuals 
on the phase transition, we use the power laws to fit the 
simulation data. Moreover, previous studies show that 
the infection rate (/3) and the virulence (a) also have a 
relation of power laws on the parameters of the infection 
rate and virulence [6l|, [62|. Yet, this relationship can be 
expressed generally with f{x) — ax^+c formula, as shown 
in Fig. [21 in which the red curve and blue curve are fit- 
ted by this function with 95% confidence intervals on the 
parameter estimates, respectively. As could be expected, 
there exist minimal and maximal critical values of the in- 
fection rate (3 on the phase transition. The open circles 
(o) and faced circles (•) indicate the simulated results 
of the minimal critical infection rate (3min and maximal 
critical infection rate Pmax, respectively. By compari- 
son between fitted curves and simulations, one can see 
there is a nonlinear trade-off, /?(t/) is a monotonically 
decreasing function of r/ and bounded by a positive con- 
stant. The dashed line in Fig. [5] show the prediction on 
the phase space from mean- field theory on TZq. From 
the classical expression of TZq , we know that the diseases 
will die out if TZq ^ 0. This has important implications, 
as it imphes that for a wide range of {f3,Tj) space there 
are two parameter domains in the {P,ti) space. One is 
domain of extinction and the other is domain of persis- 
tence. This contrast with the result obtained in spatially 
explicit model where the parameter domain is divided 
into three: domains (I) disease free; domain (II) endemic 
(see Fig. [2]). The numerical results clearly indicate (see 
Fig. [21 Inset) the validity of the prediction for the spa- 
tial thresholds to control invasion of parasites. In the 
log-log plot (Inset of Fig. [2]), the dashed fine and power 
laws (c = 0) is straight line of slope —1 shown on the 
graph. Then, the same where log(/(a;) — c) is plotted 
against log(a;). One can see that a nonzero c (nonvanish- 
ing infection rate) is better straight line. In this case, the 
nonzero c corresponds to the phase boundaries for large 
T/. The best- fit line on a log- log scale has a slope of - 
0.7065, indicating a linear dependence of log t/ on log/3. 
An understanding of this scaling behavior is lacking. 

The classical mean-field theory on the diseases spread- 
ing has shown that the number of number of secondary 
cases due to a single infected individual (referred TZq 
is maximized in the model ([2]) and therefore maximum 
transmission f3 is selected. Here, our results show that 
once the spatial structure is included the TZq is no longer 
maximized and the transmission rate is limited (see the 




FIG. 1: (Color online) Three types of spatial pattern absent 
of mobility with r/ = 0.40, Pmin ~ 0.60 and Pmax ~ 2.10. 
Here three possible states are shown: susceptible cells (white) , 
infected (red), and resistant (black), respectively, (a) Local- 
ized disease outbreaks are self-limiting in size for (3 — 0.62 
; (b) Turbulent waves for /3 = 1.20; (c) Stable spiral waves 
for (3 = 1.80. Grid size for all panels is 200 x 200. Note 
that all the patterns arise from the same initial condition. 
These spatial patterns can be produced by using supporting 
online material (see the computer code for the cellular au- 
tomata at http://7y.nuc.edu.cn/jinzhen/Englishl.aspx). The 
random number generator itself may be different if you are not 
using Linux, but it does not make a qualitative difference. 



labeled a in Fig. We note that this results are con- 
sistent with parasite evolution and extinction (6^ . In 
addition, in the spatial explicit model, the diseases can 
survive for a lower infection rate than mean-field context. 
This is a low infection rate will tend to increase the local 
density of susceptible individual around infectious. Spa- 
tial structure within cluster of individual favors a lower 
infection rate. The other way round, with higher infec- 
tion rate will tend to infect all the individuals in a cluster 
quickly. This in turn leads to relatively rapid local cluster 
extinction. 

Moreover, we also find that, for the fixed infection pe- 
riod T/, the different moving spatial patterns emerge with 



5 



infection rate j3 increasing. Figure [T] shows the typical 
snapshots of the stable spatial patterns when the param- 
eters are within domain of persistence II in Fig. [21 in 
which the results indicate that the threshold values [3min 
is equal to 0.60 and (imax equal to 2.10 for t/ — 0.40 re- 
spectively. When the infection rate is low (but just above 
the critical value Pmin), we find that the susceptible and 
infectious individual coexist, and self-organized spatial 
patterns with limited size of moving clusters. It means 
that localized disease outbreaks are self-limiting in size 
for this case (see Fig. [IJa)). With increasing infection 
rate /?, these structures grow in size and form pattern 
of moving turbulent /spirals (see Fig. [T](b) and (c)), and 
disappear for large enough /S (exceed the critical value 

Pmax ) . 
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FIG. 2: (Color online) Extinction and persistence phase tran- 
sition for (r/,/3)-parameters space. The red circles and blue 
circles are estimating from the simulation for l3max and /3min 
respectively. Red line (marked by a) is fitted by using f{x) — 
ax' + c function and a = 0.5264 (0.4456, 0.6071), 6 = -0.7324 
(-0.7647, -0.7001), and c= 1.008 (0.7124, 1.304). Blue line 
(marked by b) is fitted by using f{x) — ax'' + c function and 
a = 0.6276 (0.4384,0.8168), b = -0.6658 (-0.7287, -0.6029), 
and c = -0.3256 (-0.4643, -0.1869). Dashed line (marked 
by c) denotes the prediction obtained from mean-field theory 
on TZq. Inset shows the a nonzero c (nonvanishing infection 
rate) and zero c cases respectively on the log-log plot. 

Bartlctt observed in a series of papers [H, [g^,!!^ that 
measles in large cities had recurring outbreaks, while 
it went extinct in small communities until reintroduced 
from external sources. This means that the time to ex- 
tinction is an increasing function of the community size. 
However, the analysis by Nasell showed that Bartlett's 
approximation of the extinction time was unsatisfactory 
in an important part of parameter space (ggI . [67| . Our 
present system ([1]) is a simple individual-based stochas- 
tic SIRS model with diffusion processes [H, [sl, [tO, [zH ■ 
The goal of the analysis is to derive information about the 
quasi-stationary distribution and the time to extinction. 
Pursuing this goal leads to difficult mathematical prob- 
lems. Exact solutions cannot be found. A possible way 



to proceed is therefore to work with approximations by 
simulation. Early efforts to derive approximations of the 
expected time to extinction were shown by Nasell to lead 
to large errors (66j . In the present paper, a approximation 
of the discrete state Markov chain ([T]) is introduced. It 
led to a truncated normal distribution as an approxima- 
tion of the marginal distribution of infected individuals 
in quasi-stationarity [g^, [t^, [t^ . The resulting approx- 
imation of the expected time to extinction was rather 
coarse by a linear relation, but turned out to be an ac- 
cepted approximations for present analysis of extinction 
and persistence. 

The dependence of the average time (Tex) to extinc- 
tion on the system size of changes, N, is shown in Fig. [31 
which also demonstrates the dependence of T^x on the 
size of the populations when the parameters are within 
domain of extinction. As seen, average extinction time 
increases linearly with system size N (we test the system 
size from 10,000 to 90,000). More quantitatively, Tex is 
approximately proportional to N, which is fitted the sim- 
ulation results by using linear function, and then obtain 
the relationship as T^x = 0.055N + 758.75 are shown in 
Fig. [21 To ensure that the used averaging over just 100 in- 
dependent runs yields good statistics, we present in Fig.[3l 
the each data of Tex obtained 100 runnings for tj = 0.40 
and f3 = 0.50 (these parameters corresponding to extinc- 
tion lower domain I in Fig. [21). It should be noticed that 
previous phase transitions analysis is convincible because 
third time to extinction was used to estimate phase tran- 
sitions of persistence and extinction. One can seen that 
10 thousand time units are a long-term run for persis- 
tence on this system. Here, it is instructive to note that 
our results show that the persistence and extinction are 
independent on the population size. The opposite view 
have been obtain from the well-mixed systems in recently 
Ref. [zl] , in which they demonstrate a critical population 
size above which coexistence is likely. 
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FIG. 3: (Color online) Average time to extinction T^x versus 
system size A'^ for n = 0.40 and (3 = 0.50. (the least-squares 
linear fit of the data has a slope of 0.055 and an intercept of 
758.75 (see text for details). 
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B. Mobility promotes persistence on spatial 
epidemic model 



In fact, the individual always exhibits motion in the 
space. How does individuals mobility, in addition to non- 
linearity, affect the system's behaviors (i.e., persistence 
and extinction, spatial pattern, invasion speed and so 
on)? The insight of this important issue can be gained 
from the cellular automata model. In particular, one can 
compare the results of present paper with the previous 
works fT]. Two models are commonly used to predict 
spatial spread of a disease. The first is the distributed- 
contacts model, often described by a contact distribu- 
tion among stationary individuals. Distributed-contacts 
models are par ticularly appropriate for the study of plant 
disease [2l,lzl|> but researchers are also using the closely 
related framework of contact networks to study disease 
transmission in human population (ttI . Fzsj or animal 
species. Notice that recently developing spatial moment 
closure methods are also based on this framework jTOl.lsoj. 
The second is the distributed-infectives models, often de- 
scribed by the mobility of infected individuals. This ap- 
proach relics on the assumptions that disease is transmit- 
ted through interactions between dispersing individuals, 
and that infected individuals move in uncorrelated ran- 
dom walks. Medlock and Kot developed a distributed- 
infectives framework that uses a flexible kernel-based ap- 
proach similar to that employed in distributed-contact 
models (8]| . They found that inappropriate application 
of either the distributed-contact or distributed infectives 
approaches can generate inaccurate projections of epi- 
demic spread. Hence, here we consider a scenario that 
the transmission process involves components of both 
distributed-contacts and distributed infectives-diffusion 
of infected individuals. For instance the spread of ra- 
bies that foxes tend to be restricted to discrete home 
ranges. In Central Europe, the home range size is about 
4 km^ , but this may differ considerably between areas 
(range 2.5-16 km^) [l^, [s^]. Based on previous sub- 
section, now we investigate the individuals with a cer- 
tain form of mobility. Namely, at rate e all individu- 
als can exchange their position with a nearest neighbor. 
Each individual randomly exchanges with its eight neigh- 
bor (North, West, South, East, Northwest, Southwest, 
Southeast, and Northeast) at each time step. This is 
reasonable for the realistic system. These exchange pro- 
cesses lead to an effective mobility of the individuals. 

In this subsection we analyze the dynamics of the 
system by considering individual spatial diffusion pro- 
cess. In fact, the susceptibles' mobility will also lead 
to the spread of the disease since it changes the spatial 
distribution of the infected individuals. In present pa- 
per, we consider the mobile individuals including sus- 
ceptible, infected, and resistant. Hereafter, we refer 
the terms exchange, mobility, and diffusion to as the 
same meaning-the mobility of individuals within nearest 
neighbors. Hence, we simply refer them as "mobility". 



We refer to the mobility as following 
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where X G {S,I,R}. It has to be noted this mobility 
is not taken into account at the mean-field rate equation 
level, as well as van Ballegooijen and Boerlijst's study 
Denote L the linear size of the d-dimensional hypercubic 
lattice (i.e. the number of sites along one edge), such 
that the total number of sites reads N = L'^. Choosing 
the linear dimension of the lattice as the basic length 
unit, the macroscopic diffusion constant D of individuals 
stemming from mobility processes reads D = £d~^N~'^/'^ 
for a continuum limit [TOllssj. In this limit, a description 
of the stochastic lattice system through stochastic partial 
differential equations (SPDE) becomes feasible. But here 
we study the system's behavior in the approach of the 
cellular automata model. 

We performed cellular automata simulations of Eq. ([T]) 
with the mobility processes Eq. ([3]). The space and 
time steps were chosen the same as Subsection IIII Al In 
our simulation, the relationship between mobility rate 
and the probability of mobility apply the Gillespie al- 
gorithm H, d^. We wish to examine how the move- 
ments of individuals affect the persistence and extinction 
of epidemic in the spatially structured population. In 
particular, we will examine the role of mobility's strength 
e when the parameters are less than the critical values 
Pmm- First, we investigate the time series of the infected 
fraction. We have kept the infection period t/ fixed at 
a value r/ = 0.40, and systematically varied the mobil- 
ity rate e. In Figs. [3] and O we show typical long-time 
series of the infections for various values of the mobil- 
ity rate when the parameter (3 = 0.42 (below the critical 
values fimin ~ 0.6) and /? = 2.50 (above the critical val- 
ues Pmax ~ 2.10), respectively. Fig. UJ^a) and (b) show 
a family time series of infected fraction before and af- 
ter the mobility rate is turned on respectively, different 
e values marked by a, h, and c in Fig. SJb). We ob- 
serve in Figs, m that the mobility raises the possibility of 
shift from extinction to persistence state in the spatial 
epidemics model. Furthermore, this persistence with ir- 
regular fluctuation is around a positive equilibrium. It is 
worth noting that the large mobility rate e can promote 
the persistence of spatial epidemics in despite of the pa- 
rameter (3 below the minimal critical value Pmin- Our 
simulations indicate that the similar results hold when 
the parameter above the maximal critical value Pmax, 
but in which the persistence corresponds to low mobility 
rate e (see Fig. [SJ. 

Next, we examine the effect of mobility on the regions 
where the epidemic will extinct when the mobility rate is 
varied. By using the method developed in literature |59i] , 
we know that this scenarios can be distinguished by com- 
puting the probability Pext when the epidemic has gone 
extinct after a waiting time t oz N. For illustration, 
we have considered the flxed parameters f3 = 0.42 and 
T] — 0.40. In Fig. [HI we obtain the dependence of Pext 
on the mobility rate s. With increasing mobility rate, 
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FIG. 4: (Color online) The time-series of infected fraction 
before and after the mobility is turned on with /3 — 0.42 and 
Ti — 0.40. (a) Two independents run before mobility rate is 
turned on. (b) For different mobility rate e is turned on, line 
a: e — 10.0, line b: e — 8.0, and line c: e — 5.0. 
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FIG. 5: (Color online) The time-series of infected fraction 
before and after the mobility is turned on with /3 — 2.50 and 
Ti = 0.40. (a) Three independents run before mobility rate is 
turned on. (b) For different mobility rate e is turned on, line 
a: e — 10.0, line b: e — 8.0, and line c: e — 5.0. 



the extinction probability with a sharpened transition 
emerges at a critical value Ec = 8.0 ± 0.05 for the en- 
tire 200 X 200 lattice area explored in 10000 time-unit. 
One could see that, above Sc, the extinction probability 
Pext tends to zero as the mobility rate increase, and the 
epidemic is stable persistence (implying super-persistent 
transients ^8d\). On the other hand, below the critical 
mobility rate Ec, the extinction probability approaches 
1 for small mobility rate, and the epidemics occurs is 
unstable. One of our central results is that we have iden- 
tified a mobility threshold for persistence of epidemics. 
Without loss of generality, we also have tested the Pext 
data for the other parameters values r/ and /3, the same 
sharpened transitions have been observed from the lattice 
simulations. 

There exists a critical value Ec such that a high mo- 
bility E > Ec changes the extinction phase to persistence 
phase, while e < Ec remains extinction, leaving a uni- 
form state with only susceptible in entire space. It is 



worth noting that our above results arise from the do- 
main I of extinction phase in Fig. [5] where /3 < /3mm and 
absence of mobility. Therefore, our results also indicate 
that the stronger mobility can increase the parameter re- 
gion where the epidemic persistence occurs. This may 
be completely explained from the local and global infec- 
tion [6l|] . The higher mobility leads to a increasing of the 
global contact among individuals in the lattices. 




mobility, e 



FIG. 6: (Color online) The critical mobility Ec with different 
system sizes. Mobility above the value Sc induces persistence; 
while it is extinction below that threshold. The data obtained 
from lattice simulations of the system after long temporal 
development (that is, at time t oc N) and for different values of 
e with P = 0.42 and tj = 0.40. Here, we have considered the 
extinction probability Pext starting with randomly distributed 
individuals on a 200 x 200 square lattices and 100 x 100 square 
lattices, respectively. 

Reichenbach et al had assessed the effects of mobility 
on the spatial pattern of rock-paper-scissors games |59j. 
Results demonstrated that a critical influence of mobil- 
ity on species diversity: when mobility exceeds a certain 
value, biodiversity is jeopardized and lost. In contrast, 
below this critical threshold all subpopulations coexist 
and an entanglement of traveling spiral waves forms in 
the course of time. In our study, we incorporate the sce- 
narios outlined in Refs. 59, 87] and also investigate a spa- 
tial pattern induced by mobility when the parameters are 
within domain II of Fig. [2l We perform extensive com- 
puter simulations of the stochastic system (see Eqs. ^ 
and ([3])) and typical snapshots of the steady states are 
shown in Fig. [Tj With increasing mobility e, these mov- 
ing spiral structure grow in size and saturate at a certain 
value, but they don't disappear with respect time. This is 
different from rock-paper-scissors games where the spiral 
wave will disappear for large mobility [59[. As shown in 
Fig. [3 the spirals' wavelength A rises with the individu- 
als' mobility increased. When the parameters are within 
domain I, we have not observed the spatial pattern emer- 
gence for varied mobility rate e. However, we find that 
the high mobility lead to entire disease outbreaks instead 
of localized disease outbreaks (see Fig. |8]). 

From the above observations we know that the mobil- 



ity of individuals can dramatically affect the epidemics 
extinction/persistence and spatial pattern. We need to 
point out here that, although our results in this study 
come from cellular automata's simulation, it is also useful 
for assessing the dynamics of deterministic systems, such 
as stochastic PDEs, spatial correlation models, ordi nary 
difference equation, and so forth [see [H, [sl, [TO, [H, [s^ . 
By modeling the interaction among individuals, we are 
able to understand the role of spatial mixed (cause by 
mobility) in invasion dynamics without the need for com- 
plex mathematical methods. Whether an organism can 
successfully invade and persist in the long-term is depen- 
dent on many factors. Here, we have identified a mobility 
factor for epidemic persistence and invasion, which is dis- 
tinguished from the nonlinear dynamics of the epidemic 
models. 



(a) 




FIG. 8: (Color online) Spatial pattern in the lattice simulation 
for mobility above the value Sc induces persistence with /3 = 
0.42 and tj = 0.40. Here, we have considered the spatial 
pattern starting with randomly distributed individuals on a 
square lattice, (a) e = 8.0; (b) e = 12.0. 



scales. To quantify the traveling behavior of individuals, 
we consider the linked dynamics of two host communities 
each of size N, where individuals possess mobility in each 
patch. Spatial coupling is the critical parameter deter- 
mining phase coherence and spatial synchrony Jl8| and 
the persistence of host-pathogen systems [H, l93j |. For 
simplicity, we assume that N is constant through time 
for each patch. The typical schematic for two patches is 
shown in Fig. [S] 



FIG. 7; (Color online) The effect of mobility on spiralling 
patterns. We show snapshots obtained from lattice simula- 
tion of typical states of the system after long temporal devel- 
opment and for different values of e starting with randomly 
distributed individuals on a square lattice, (a) e — 0.01; (b) 
e = 8.05. Parameters are tj = 0.40 and (3 — 1.80. 



Spatio-temporal synchronous dynamics of 
spatial epidemics 




patch J 



o 



As individuals travel around the world for human, as 
well as migrate for animals, the disease may spread from 
one place to another. The spatial spread of epidemics has 
been much studied, particularly with respect to pandemic 
invasion waves [1, |a, |^, [H, Ho, HH ■ Simulation models 
incorporating transportation have generated important 
insights into the spread of epidemics. However, the key 
underlying relationship between human movement and 
disease spread has not been verified across wide spatial 



FIG. 9: (Color online) Typical spatial structured model of 
the two-patches spread of epidemics (schematic). The two 
patches i and j are connected by migration processes. Each 
patch contains a population of individuals who are character- 
ized with respect to their stage of the disease (e.g. susceptible, 
infected, resistant), and identified with a different color in the 
picture. 

Consider now, two patches with sizes 200 x 200 cells 
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each other, each representing a community or suburb, 
and suppose a smaU number of infective individuals in 
each community at initial time. The individual can pass 
from one community to the other. This migration or the 
individual traveling is randomly chosen for some given in- 
dividuals moving from a community to another commu- 
nity in the large metapopulation systems. The strength 
of coupling between the two communities is measured by 

LR — no. of migration/size of community (4) 

with respect to unit time. Although the approach is 
only a rough approximation of the dynamic nature of 
transportation between two communities, it neverthe- 
less capture some of the essential elements [1, |93| • Fig- 
ure [10] shows a simulation results with different coupling 
strength between two communities, in which the mobil- 
ity within patch is ignored. In this case, the infectives 
in each community are seen to irregular oscillate and lo- 
cal anti-phase for weak coupling strength (see Fig. [TUTa') 
and (b)). Upon increasing the coupling strength between 
the communities, the suburbs show in-phase synchroniza- 
tion jl06l | (see Fig. [TUf c)). Notice that the similar results 
have been obtained by using SIRS network model [1, [2l[ . 
As previous subsection's results show that, at the sin- 
gle patch level, the epidemic behavior on the global [l07j | 
scale is also determined by the mobility of individuals. 
In particular, the effects due to the finite size of com- 
munities and the stochastic nature of the mobility might 
have a crucial role in the problem of resurgen t ep idemics, 
extinction and eradication [also see[9^[95l. l96l[97| . There- 
fore it is important to consider the effect of mobility rate 
on synchronization for the metapopulation system. In 
Fig. HH we report the effect of different mobility rate on 
the time-series of infected fraction by the lattice simula- 
tions. Figure. [TT] provides a clear evidence of the spatio- 
temporal synchrony in metapopulation system indepen- 
dent of the self-mobility, but it depends on the coupling 
strength among patches. One can draw a conclusion that 
mobility within each patch does not play a signification 
role, and spatial structure within patches can be ignored 
when investigate the coupling patches. 

One of our main goals is to study the manner in which 
two or more communities or suburbs synchronize in time, 
as recurring waves of infection sweep through their re- 
spective populations by means of cellular automata sim- 
ulation. This is a fascinating outcome of spatial dynam- 
ics whereby a few migrating infective individuals have 
the potential to spread the epidemic from one commu- 
nity to the other, giving rise to synchrony in the long 
term. For instance, the time-series of measles infections 
in various cities of England (1944-1958) are reported in 
literatures [1, [3 ■ The Birmingham and Newcastle 
appear to synchronize in-phase together while Cambridge 
and Norwich are clearly out of phase by 180 degree (see 
Fig. [12] in Appendix and Fig. 2 in literature [ll, \9^). 
Here, we simulate a simple coupled spatial community 
models to gain insights about the coupling strength. 
These results provide useful insights for the basic the- 
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FIG. 10: (Color online) The time series of infected fraction 
in two coupled communities with e = 0, /? = 1.80, and 
Ti = 0.40. (a) With weak coupling, LR = 0.001, rapidly 
anti-phase synchronization predominates; (b) With middle 
coupling, LR = 0.01, remaining anti-phase synchronization 
predominates; (c) With strong coupling, LR = 0.1, epidemic 
peaks, and the troughs between them, rapidly synchronize 
in-phase. 
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FIG. 11: (Color online) The time series of infected fraction 
in two coupled communities with different mobility rate in- 
side each patch, (a) Low mobility rate, e — 0.5, anti-phase 
synchronization predominates; (b) With hight mobility rate, 
£ — 4.5, remaining anti-phase synchronization predominates. 
The other parameters are same as in Fig. llOf b') . 



oretical understanding of mechanistic epidemic models 
in complex multi-patches environments, which can then 
be used to build more realistic data-driven large-scale 
computational approaches for real case scenarios and 
spatially targeted control measures. For instance, the 
spread of influenza shows spatial synchrony in 49 states 
of United Sates 3 ■ 
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As a summary, from the simulations of individuals' 
spreading processes among lattices with randomly dis- 
tributed in the space, we can identify different effects of 
parameters describing the spatial structure and the tem- 
poral developments of disease for individuals, such as r/, 
/3, e, LR and N, on the dynamical behavior-persistence, 
extinction and synchrony. We change every parameter 
which we consider can affect the spreading process. Al- 
though the short-term behavior (persistence and extinc- 
tion) may depend on system size N and details of spa- 
tial distribution, the long-term behavior mainly depends 
on three parameters: the infection period t/, infection 
rate /3, and mobility e for single community. The two 
former parameters determine what types spatial pattern 
emerge, and the last dramatically increase the parame- 
ters region where the epidemics pandemic occurs in the 
space. Moreover, we investigate the coupling between 
multi-community, in which the two former parameter t/ 
and f3 reflect the nature of the developing and pandemic 
period of disease for individuals, the coupling strength 
LR is a parameter related to the population density, the 
frequencies of people contacts, and the extent of people 
traveling between different places or cities [9^. These 
results may be helpful for the analysis of the spreading 
processes of diseases in space. 

The possibility of spiral waves, self-organized spatial 
pattern, spotted pattern, traveling wave as well as trade- 
offs in spatial epidemics and host-parasitoid b y lo cal dis- 
persal abilit ies has long been recognized P, [l,Ill, [H, [H, 
[23, I10G| | . What is new here is that the trade-off be- 
tween persistence and extinction have been understood 
in the spatial structure epidemics model. The discrete 
character of the individuals involved in the reactions ([T]) 
and the mobility processes ([3]) are responsible for intrinsic 
stochasticity arising in the system. Individuals' mobility 
as well as intrinsic noise have crucial influence on the self- 
formation of spatial patterns. The analytical expressions 
for the spirals' wavelength as function of mobility can 
be determined by means of a complex Ginzburg-Landau 
equation (CGLE) obtained by recasting the PDE derived 
from interacting particle approach [53 . pTOj . Hence, to 
qualitatively explain these findings, a profound under- 
standing is still desirable and could motivate further in- 
vestigations. 



APPENDIX A 

We give here the time series of weekly measles case 
reports for Birmingham and Newcastle between 1948 and 
1969, which exhibits in-phase pattern. 
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FIG. 12: (Color online) The time series of weekly measles 
case reports for Birmingham and Newcastle between 1948 
and 1969. Data made available by Professor Benjamin Bolker 
(www.zoo.ufl.edu/boIker/measdata.html). 
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